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Abstract : In the present study, autoregressive integrated moving average (ARIMA) methodology has been applied for modeling 
and forecasting of yearly area and production of rice in India. Rice production data for the period of 1950-195 1 to 2014-2015 of 
India were analyzed by time-series methods. Autocorrelation and partial autocorrelation functions have been estimated, which 
have led to the identification and construction of ARIMA models, suitable in explaining the time series and forecasting the future 
area and production.The diagnostic checking has shown that ARIMA (1,0,1) and ARIMA (0, 1 , 1 ) is appropriate for rice area and 
production. The forecasts from 2015-2016 to 2024-2025 were calculated based on the selected model. The forecasting power of 
autoregressive integrated moving average model was used to forecast rice area and production for ten leading years. This 
projection is important as it helps to inform good policies with respect to relative production, price structure as well as consumption 
of rice in the country. 
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Introduction 

Paddy is the most important and extensively grown 
food crop in the world. It is the staple food of more than 
60 per cent of the world population. It is of special 
importance for the nutrition of large reaches of the 
population in Asia. As a result, it plays a pivotal role for 
the food security of over half the world population. It is 
also a central component of the culture of a number of 
communities. For those reasons, rice is considered as a 
“stratregic” commodity in many countries, both 
developed and developing, and has consequently remained 
subject to a wide range of government controls and 
interventions. India has the largest area under paddy in 


the world and ranks second in the production after China. 
Country has also emerged as a major rice consumer. 

Rice is cultivated in more than 100 countries in the 
world. In 2015-16, China produced some 145.5 million 
metric tons of milled rice, making the country the world’s 
largest producer of rice followed by India (103.5 million 
tons). Together these countries, accounted about half of 
world rice area and production. Indonesia (36.3 million 
tons), Bangladesh (34.6 million tons), Vietnam (28.2 
million tons), Thailand (16.4 million tons), Burma (12.2 
million tons), Philippines (11.5 million tons), Brazil (8 
million tons), and Japan (7.9 million tons) are the other 
major rice producing countries. ( www.statista.com ). In 
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case of productivity, Egypt ranks first with 9086 kg/ha 
followed by USA (7037 kg/ha), Japan (6702 kg/ha) and 
Korea Rep (6592 kg/ha). According to 4 th advance 
estimate, India accounted for 439.49 lakh hectares area 
with production level of 106.54 million tons of rice. The 
productivity of rice in India is 2424 kg/ha which is less 
than other rice producing countries. 

Forecasts have traditionally been made using 
structural econometric models. Concentration have been 
given on the univariate time series models known as auto 
regressing integrated moving average (ARIMA) models, 
which are primarily due to world of Box and Jenkins 
(1970). These models have been extensively used in 
practice for forecasting economic time series, inventory 
and sales modeling (Brown, 1959 and Holt et al, 1960) 
and are generalization of the exponentially weighted 
moving average process. Several methods for identifying 
special cases of ARIMA models have been suggested 
by Box-Jenkins and others. Makridakis et al. (1982) and 
Meese and Geweke (1982) have discussed the methods 
of identifying univariate models. Among others Jenkins 
and Watts (1968); Yule (1926 and 1927); Bartlett (1964); 
Quenouille (1949); Ljune and Box (1978) and Prindycke 
and Rubinfeld (1981) have also emphasized the use of 
ARIMA models. 

In this study, these models were applied to forecast 
the area and production of rice crop in India. This would 
enable to predict expected rice area and production for 
the years from 2015 onwards. Such an exercise would 
enable the policy makers to foresee ahead of time the 
future requirements for grain storage, import and/or 
export of rice thereby enabling them to take appropriate 
measures in this regard. 

Material and Methods 

Time series data were used for the study. The data 
were obtained from the website of Directorate of 
Economics and Statistics, Department of Agriculture and 
Cooperation, Ministry of Agriculture from 1950-51 to 
2014-15. Box and Jenkin (1976) linear time series model 
was applied on the data. This model is commonly known 
as autoregressive integrated moving average model 
(ARIMA Model). 

One of time series models which is popular and 
mostly used in ARIMA model. ARIMA (p, d, q) model is 
a mixture of autoregressive (AR) model which shows 
that there is a relation of a value in the present (Z) and 
value in the past (Z t k ), added by random value and 


moving average (MA) model which shows that there is 
a relation between a value in the present (Z t ) and residuals 
in the past 

(Z, k k = 1,2...) 

with a non-stationary data pattern and d differencing 
order. The form of ARIMA (p,d,q) is: 

® P (B) (1- B ) d Z t = Q q (B)a t 

where, p is AR model order, q is MA model order, d 
is differencing order and: 

<t> (B) = (1- <t> B - <t> £ 2 - BP) 

p v v 1 2 P 

3> (B) = (1- <J>, B - 4> B 2 - ...-<P B « ) 

v l 2 q 

Generalization of ARIMA model for a seasonal 
patter data, which is written as: 

ARIM ( p , d , q){ P , D , Qf= <£, (B) <£ (B 2 ) (1 -B ) d (1-B ) (1 -B° ) D 
Z t =e (B) & Q (B s )at 

where, s is seasonal period. 

(B s ) = (l-c^B -<t> 7 B 2 - B p ) 

e^BMi-e, B-e 2 b 2 -...-qb«) ' 

Model identification : 

To determine whether the series isstationary or not 
we considered the graph of ACF. If a graph of ACF 
cuts of fairly quickly or dies down fairly quickly, then the 
time series value should be considered stationary. Model 
for non-seasonal series are called autoregressive 
integrated moving average model, denoted by ARIMA 
(p,d,q). Here p indicates the order of the autoregressive 
part, d indicates the amount average of difference and q 
indicates the order of the moving average part. If the 
original series is stationary, d = 0 and the ARIMA models 
reduce to the ARMA models. 

The difference linear operator (A), denoted by: 

AY t = Y t - Y ,i = Y t - BY t = ( 1 - B )Y t 

The stationary series: 

W= A d Y=(l-B ) d Y= p +0 ? ( B)e t 

or <P p ( B )W= p +0/ B)z t 

Model estimation and checking : 

Estimate the parameters for a tentative model has 
been selected. The derived model must be checked for 
adequacy by considering the properties of the residuals 
whether the residuals from an ARIMA model is normal 
and randomly distribution. An overall check of the model 
adequacy is provided by Ljung-Box Q statistics. The test 
statistics Q is given in equation below: 

n Tk (e) 2 

Qm= n (n + 2)I n k m-r 

k-1 
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where, r k (e) = The residual autocorrelation at lag K. 
n = The number of residuals 
m = The number of time lags includes in the test. 
If the p-value associated with the Q statistics is 
small (p-value <“), the model is considered inadequate. 
The analysts should consider a new or modified model 
and continue the analysis until a satisfactory model has 
been determined. 

Results and Discussion 

The maximum cultivation Area of rice was 45.16 
million hectares in 1999-00 and was minimum 29.83 
million hectares in 1951-52. For production, the maximum 
production of rice in India was obtained in 2013-2014 
year (106.65 million tons) and minimum in 1950-1951 
year (20.58 million tons). 

Last 65 years data of cultivation area and production 


of rice in India was used for modeling purpose. In model 
specification, we looked at the plots of auto-correlation 
function (ACF) for rice cultivation areas (Fig. 1) and 
production figures (Fig. 2). Also, partial auto correlation 
function (PACF) for rice cultivation area (Fig. 3) and 
production (Fig. 4). Auto correlation function indicated 
the order of the auto regression compounds “q” of the 
model while the partial correlation function gave an 
indication for the parameter p. The ACF and PACF of 
the residuals (Fig. 5 and 6) also indicate ‘good fit’ of the 
model. 

The time series plot (Fig. 7 and 8) of cultivation 
areas and production showed an increasing trend. ACF 
of both series showed non-stationary as ACF did not fall 
as quickly as the log K increased. To check the further 
stationary, second difference of the original series for 
cultivation area and production was taken. The auto 
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Fig. 1 : Autocorrelation function of rice cultivated area 
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Fig. 2 : Autocorrelation function of rice production 
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Fig. 3 : Partial autocorrelation function for rice cultivated area 
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Fig. 4 : Partial autocorrelation function for rice production 
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Fig. 5 : ACF of residuals of fitted ARIMA model 
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Residual ACF Residual PACF 
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Fig. 6 : PACF of residuals of fitted ARIMA model 
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Fig. 7 : Time series plot for cultivated area 
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Fig. 8 : Time series plot of production area 


Internat. J. agric. Sci. 1 June, 2017 I Vol. 13 I Issue 2 1 208-214 



Hind Agricultural Research and Training Institute 



PREMA BORKAR 


correlation formulation of second series and correlogram 
shows some more stationary than that of the first 
different. The corellogram of the auto correlation function 
of first difference series showed that the auto correlation 
function falls finally after lag 1 for area and production, 
hence the respective values of the parameter “q” decided 
to be 1 . 

PAC function of the first differenced series of 
the cultivation area and production was used to determine 
parameter “p”. Thus, we chose “p” to be 1 and 0 for 
cultivation area and production respectively which gave 
good results consequently, the respective value of p,d,q 
were determined for ARIMA, that is ARIMA (1,0,1) 
and ARIMA (0,1,1). 

Model estimation : 

ARIMA (1,0,1) and (0,1,1) model were estimated 
using SPSS and estimation of the models for the rice 
area and production are given in Table 1 and 2. Goodness 
of fit of the model given in as the diagnostic check of the 
estimated model. 

Residual analysis : 

The time series plot of the residual cultivated area 
and production data showed scattered trend, therefore, 


models were fitted properly by residual analysis. 

For normality test, Shapiro-wilk test was used. The 
test was significant and assumption of normality was 
accepted. Since the series fitted shows normality, the 
model is a good fit. 

Forecast of area and production : 

ARIMA (1,0,1) and ARIMA(0,1,1) were taken for 
10 years ahead and forecasts for rice area and production 
are given in Table 3 at 95 per cent confidence interval 
values. 

For 2015-16, a forecast of rice area was 42.82 
million ha with lower and upper limit of 40.27 and 45.37 
million ha, respectively. A rice area forecast for the year 
2025 was 42.45 million ha with lower and upper limit of 
36.82 and 48.08 million ha, respectively. Forecast of rice 
production showed an increasing trend. For 2015-16, a 
rice production was about 1 1 1 .23 million tons with lower 
and upper limit 89.4 and 133.06 million tons. The rice 
production forecast for the year 2025 is about 139.25 
million tons. 

Conclusion : 

In our study, the developed model for rice area and 
production was found to be ARIMA (1,0,1) and (0,1,1). 


Table 1: Estimate of area parameters 

Type 

Co-efficients 

S.D. 

t-ratio 

Constant 

3.601 

0.185 

19.452 

AR1 

0.993 

0.015 

66.602 

MAI 

0.331 

0.127 

2.605 


Table 2: Estimate of production parameters 

Type 

Co-efficients 

S.D. 

t-ratio 

Constant 

0.025 

0.003 

7.657 

AR1 

0.998 

0.001 

9.654 

MAI 

0.732 

0.092 

7.969 


Table 3 : Forecasts for rice area and production (2015-16 to 2024-2025) 

Years 

Area (Million hectares) 

Production (Million tonnes) 

2015-2016 

42.82+2.55 

111.23+21.83 

2016-2017 

42.78±3.06 

1 14.04±23.23 

2017-2018 

42.74±3.50 

116.92±24.66 

2018-2019 

42.69±3.89 

119.88±26.12 

2019-2020 

42.65+4.24 

122.91+27.62 

2020-2021 

42.61+4.56 

126.02+29.16 

2021-2022 

42.57±4.85 

129.20±30.74 

2022-2023 

42.53±5.13 

132.47±32.36 

2023-2024 

42.49±5.38 

135.82±34.03 

2024-2025 

42.45+5.63 

139.25+35.75 
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The forecasts of rice area and production, lower control 
limits (LCL) and upper control limits (UCL) are 
presented in Table 3. The validity of the forecasted values 
can be checked when the data for the lead periods 
become available. ARIMA model being stochastic in 
nature, it could be successfully used for modeling as well 
as forecasting the rice area and production of India. The 
model demonstrated a good performance in terms of 
explaining variability and predicting power. The supply 
projection of an agricultural commodity especially rice 
plays a vital role in the adjustment of supply to demand. 
These projections help the government to make policies 
with regards to relative price structure, production and 
consumption patterns and also, to establish relationship 
with other countries of the world. 
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